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There are now two dominant models of how stars form: gravitational collapse 
theory holds that star-forming molecular clumps, typically hundreds to thou- 
sands of M© in mass, fragment into gaseous cores that subsequently collapse 
to make individual stars or small multiple systems^"^. In contrast, competitive 
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Q ■ accretion theory suggests that at birth all stars are much smaller than the typi- 
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^ I cal stellar mass (~ 0.5 Mq), and that final stellar masses are determined by the 



subsequent accretion of unbound gas from the clump . Competitive accretion 



5t , models explain brown dwarfs and free-floating planets as protostars ejected from 
star-forming clumps before accreting much mass, predicting that they should 
lack disks, have high velocity dispersions, and form more frequently in denser 
clumps^"^^. They also predict that mean stellar mass should vary within the 
Galaxy*. Here we derive a simple estimate for the rate of competitive accretion 
as a function of the star-forming environment, based partly on simulations^^, 
and determine in what types of environments competitive accretion can occur. 

1 



We show that no observed steir-forming region produces significant competitive 
accretion, and that simulations that show competitive accretion do so because 
their properties differ from those determined by observation. Our result shows 
that steirs form by gravitational collapse, and explains why observations have 
failed to confirm predictions of the competitive accretion scenario. 

In both theories, a star initially forms when a gravitationally bound gas core collapses. 
The crucial distinction between them is their prediction for what happens subsequently. In 
gravitational collapse, after a protostar has consumed or expelled all the gas in its initial 
core, it may continue accreting from its parent clump. However, it will not accrete enough 
to substantially change its mass^^'^^. In contrast, competitive accretion requires that the 
amount accreted after the initial core is consumed be substantially larger than the proto- 
stellar mass. We define /m = m^tdyn/''^* as the fractional change in mass that a protostar 
of mass m* undergoes each dynamical time idyn of its parent clump, starting after the ini- 
tial core has been consumed. Gravitational collapse holds that /m <^ 1, while competitive 
accretion requires /m ^ 1- 

Consider a protostar embedded in a molecular clump of mass M and mass-weighted 
one-dimensional velocity dispersion a. Competitive accretion theories usually begin with 
seed protostars of mass m* 0.1 Mq^"''', so we adopt this as a typical value. We consider 
two possible geometries: spherical clumps of radius R and filaments of radius R and length 
L, L ^ R. These extremes bracket real star-forming clumps, which have a range of aspect 
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ratios. The virial mass for (spherical, filamentary) clumps is 

and the virial parameter is a^ir = Mvir/M^^'^^. The dynamical time is tdyn = R/o'. 

First suppose that the gas that the protostar is accreting is not collected into bound 
structures on scales smaller than the entire clump. Since the gas is unbound, we may neglect 
its self-gravity and treat this as a problem of a non-self-gravitating gas accreting onto a 
point particle. This process is Bondi-Hoyle accretion in a turbulent medium, which gives an 
accretion rate^^ 

_(Gm )2 

where p is the mean density in the clump. The factor 0bh represents the effects of turbulence, 
which we estimate in terms of a, and R in the Supplementary Information.^^ From Q 
and the definitions of the virial parameter and the dynamical time, we find that accretion 
of unbound gas gives 

/m-BH = f 14.4, 3.08^) (/.BHtt^i? f (3) 
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for a (spherical, filamentary) star- forming region. From this result, we can immediately see 
that competitive accretion is most effective in low mass clumps with virial parameters much 
smaller than unity. 

Table |21 shows a broad sample of observed star-forming regions. None of them have a 

value of /m-BH near unity, inconsistent with competitive accretion and in agreement with 

gravitational collapse. Note that the Bondi-Hoyle rate is an upper limit on the accretion. 

If stars are sufficiently close-packed, their tidal radii will be smaller than their Bondi-Hoyle 
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radii, and the accretion rate will be lower^. Also note that radiation pressure will halt 
Bondi-Hoyle accretion onto stars larger than ~ 10 Mq^^. 



The second possible way that a star could gain mass is by capturing and accreting other 
gravitationally bound cores. We can analyze this process by some simple approximations. 
First, when a star collides with a core it begins accreting gas from it, causing a drag force^*. 
If drag dissipates enough energy, the two become bound. We can therefore compute a critical 
velocity below which any collision will lead to a capture and above which it will not. Second, 
cores and stars should inherit the velocity dispersion of the gas from which they form, so we 
assume they have Maxwellian velocity distributions with dispersion a. The true functional 
form may be different, but this will only affect our estimates by a factor of order unity. 
Third, we neglect the range of core sizes, and assume that all cores have a generic radius i?co 
and mass M^o- Competitive accretion requires M^o < m*, so we take M^o — m*, which gives 
the highest possible capture rate. Finally, we make use of an important observational result: 
cores within a molecular clump have roughly the same surface density as the clump itself 
S = M{ttR'^ ,2RL)~^ for (spherical, filamentary) clumps. This enables us to compute the 
escape velocity from the surface of a core in terms of the properties of the clump. 



With these approximations, it is straightforward to compute the amount of mass that 
a protostar can expect to gain by capturing other cores. In the Supplementary Information, 
we show it is 





m— cap 



(0.42, 0.36) [4 + 2ui, - (4 + 7.32 ul,) exp (-1.33 m^) 



(5) 
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where 0co is the fraction of the parent clump mass that is in bound cores and Uesc = Vcsc/o'. 
Surveys generally find core mass fractions of 0co ~ 0.1^°"^^, so we adopt this as a typical 
value, giving the numerical values of /m-cap shown in Table |21 As with /m-BH, all the 
estimated values are well below unity. 

If we let /in = /m-BH + /m-cap, then we can use our simple models to determine where 
in parameter space a star- forming clump must fall to have /m > 1- For simplicity, consider a 
spherical clump with fixed 0bh = 5 and (pco = 0.1 (typical values for observed regions), and a 
seed protostar of mass = 0.1. In this case, both /m-BH and /m-cap are functions of a^^^M 
alone, and we find /m > 1 for al^j.M < 8.4 Mq. The functional dependence is more complex 
if we include filamentary regions and allow 0bh and 0co to vary, but the qualitative result is 
unchanged. Observed star- forming regions have ctvir ~ 1 and M ^ 10^ — 10^ Mq"^^, which 
produces /m <^ 1- No known star-forming region has al-^j.M small enough for competitive 
accretion to work. Thus, the cores from which stars form must contain all the mass they 
will ever have, which is the gravitational collapse model. 

Our simple estimate of /m is consistent with simulations of competitive accretion as 
well, and explains why competitive accretion works in the simulations. All competitive 
accretion simulations have virial parameters ctvir -C 1. In some cases the simulations start in 
this condition^' with Ovir ~ 0.01 as a typical choice. In other cases, the virial parameter 
is initially of order unity, but as turbulence decays in the simulation it decreases to ^ 1 in 
roughly a crossing time'^'^'^°'^^. Once competitive accretion gets going in these simulations, 
they have reached a^i,- <^ 1 as well. In addition, many of the simulations consider star- 
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forming clumps of masses considerably smaller than the ~ 5000 Mq typical of most galactic 
star formation^^, with M ^ 100 Mq not uncommon. Consequently, the simulations have 
a^i^M ^10 Mq, which explains why they find competitive accretion to be important. Note 
that simulations where turbulence decays will have 0bh ~ 1, rather than the typical value of 
0BH = 5 we have used for real regions, but this does not substantially modify our conclusions. 

Three other aspects of the simulations even further increase their estimate of /m- First, 
the Bondi-Hoyle radius of a 0.1 Mq seed protostar in a typical clump is only 5 AU, a smaller 
scale than any of the competitive accretion simulations resolve. This under-resolution may 
enhance accretion^^. Second, small virial parameters lead most of the mass to collapse to 
stars, giving 0co ~ 0.5 — 1 after a dynamical time, and also tend to make the cloud fragment 
into smaller pieces, lowering M. Third, rapid collapse leaves no time for large cores to 
assemble. For example, one simulation of a ~ 1000 Mq clump produces no cores larger than 
1 Mq^, inconsistent with observations that find numerous cores more massive than this in 
similar regions^^'^^. With no large cores, large stars can only form via competitive accretion. 

Thus, our results are consistent with the simulations, but they show that the simula- 
tions are not modelling realistic star-forming clumps. One might argue that all clumps do 
enter a phase with ccyir <S 1 that occurs rapidly and has therefore never been observed, but 
that most stars are formed during this collapse phase. In this scenario, though, protostars 
associated with observed star-forming regions should have systematically lower masses than 
the field star population, since they were formed before the collapse phase when competitive 
accretion might work. One would also expect to see a systematic variation in mean stellar 
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mass with age in young clusters, corresponding to cluster evolution into a state more and 
more favorable to competitive accretion. This is not observed. 

We hypothesize that the primary problem with the simulations, the reason they evolve 
to CKvir <^ 1, is that they omit feedback from star formation. Recent observations of pro- 
tostellar outflow cavities show that outflows inject enough energy to sustain the turbulence 
and prevent the virial parameter from declining to values much less than unity^^. Another 
possible problem in the simulations is that they simulate isolated clumps containing too 
little material. Real clumps are embedded in molecular clouds, and large-scale turbulent 
motions in the clouds may cascade down to the clump scale and prevent the turbulence from 
decaying. A third possibility is that turbulence decays too quickly in the simulations because 
they do not include magnetic fields and their initial velocity fields, unlike in real clumps, are 
balanced rather than imbalanced between left- and right- propogating modes^^. 

One implication of our work is that brown dwarfs need not have been ejected from their 
natal clump, so their velocity dispersions should be at most slightly greater than those of 
stars, and their frequency need not change as a function of clump density. This also removes 
a discrepancy between observations showing that brown dwarfs have disks^^ and theoretical 
models of their origins. Another conclusion is that the mean stellar mass need not vary from 
one star-forming region to another as competitive accretion predicts, removing a discrepancy 
between theory^ and observations that have thus far failed to find any substantial variation in 
typical stellar mass with star-forming environment. In the gravitational collapse scenario, the 
mean stellar mass may be roughly constant in the Galaxy, but may vary with the background 
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radiation field in starburst regions and in the early universe^. 
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Name 


Mass 


Type 


M (Mq) 


R (pc) 


L (pc) 


/I 1 \ 

cr (km s J 


L1495 P° 


Low 


Sph. 


410 


2.1 




0.58 


L1495 Ipo 


Low 


Sph. 


950 


2.4 




0.67 


11709^^ 


Low 


FiL 


140 


0.23 


3.6 


0.48 


L1755^^ 


Low 


FiL 


171 


0.15 


6.3 


0.53 


W4423 


High 


Sph. 


16000 


0.35 




3.9 


W75(OH)23 


High 


Sph. 


5600 


0.25 




3.5 




High 


FiL 


5000 


0.25 


13 


1.41 


N-fil^^ 


High 


Fil. 


16000 


2.3 


88 


1.54 



Table 1: Sample star-forming regions. Sph. = sphericaL FiL = Filamentary, /-fil = Orion 
integral filament, N-fil = Orion North filament. For L1495 I and II, the data are from Kramer &i 
Winnewisser's ^^CO observations, and the masses are Kramer &i Winnewisser's Mqo- For W44 
and W75(OH), the data are from Plume et al.'s CS J = 5 — > 4 observations, and the masses are 
Plume et al.'s M^. Note that, since CS J = 5 — > 4 is a very high density tracer, it biases the 
results to small virial parameters by excluding low density parts of the clump. 
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Name ctyir 4>bh fm — BH ^esc /m— cap 

L1495 I 2.0 2.4 0.0022 0.28 0.0015 

L1495 II 1.3 3.0 0.0026 0.28 0.0015 

L1709 2.8 0.93 0.0042 0.44 0.0072 

L1755 4.8 0.54 0.0017 0.40 0.0052 

W44 0.39 6.4 0.0038 0.25 0.0010 

W75(OH) 0.63 5.2 0.0034 0.26 0.0011 

/-fil 2.4 4.1 0.0023 0.26 0.00097 

N-fil 6.2 5.7 0.0001 0.11 0.00003 
Table 2: Computed properties of sample star-forming regions, /-fil = Orion integral 
filament, N-fil = Orion North filament. 
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